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Abstract 

The study of possible new physics signals in global event properties in 
pp collisions in the TeV energy domain is extended from full phase- 
space to rapidity intervals experimentally accessible at LHC The elbow 
structure in the total multiplicity distribution predicted in full phase- 
space is clearly present also in restricted rapidity intervals, leading to 
very strong charged particle correlations. It is also found that energy 
densities comparable to those reached in heavy ion collisions at RHIC 
could be attained in pp collisions at LHC. 
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Introduction 

The search for signals of new physics in global event properties in pp collisions 
at LHC in the framework of the weighted superposition mechanism of different 
classes of minimum bias events as components led us to explore the possibility of 
the existence at 14 TeV cm. energy of a new class of events []]]. The aim of the 
present paper is to extend our study from full phase-space (FPS) to pseudo-rapidity 
intervals which will become experimentally accessible with the LHC detectors. 

The paper is organised as follows. In Section 1 are summarised results of 
on general properties of soft and semi-hard components in pseudo-rapidity intervals 
in extrapolated scenarios based on the knowledge of the GeV energy domain. In 
Section 2 the existence at 14 TeV cm. energy in pp collisions of a third class of events 
in addition to the soft and semi-hard ones is postulated following our findings on 
the same class in FPS and its general properties in pseudo-rapidity intervals are 
discussed. For completeness, a comparison with Pythia Monte Carlo calculation 
results is performed where appropriate. 



1 Antecedents: global properties in full phase-space 

In jif the total n-charged particle multiplicity distribution in pseudo-rapidity inter- 
vals \r)\ < T] c as a function of cm. energy \fs was written as follows: 

PniVc \/s) = oi soit (\/s)P^ B (n soft (r] c , y/s), k solt (r] c , \/s))+ 

(1 - « so ft(v^))^„ PaNB (nsh(^, v^), k Bh (rj c , Vi)), (1) 

where P BaNB (n, k) are Pascal (negative binomial) multiplicity distributions with char- 
acteristic parameters n, the average charged multiplicity, and k, related to the 
dispersion D by k^ 1 = (D 2 — n)/n 2 , in the pseudo-rapidity interval \r}\ < r) c at 
cm. energy y/s for the soft and semi-hard components, and a so ft is the fraction of 
soft events with respect to the total number of events. The most appealing result 
has been the decrease of the average number of clans of the semi-hard component 
N s h(Vc) at fixed r\ c as y/s increases (already seen in FPS both in the strong-KNO- 
scaling- violating scenario (where fe^Li — l n s ) an d in the QCD-inspired one (where 
A;^ 1 ~ 1/Vln s), together with the corresponding increase of the average number of 
particles per clan. In addition, the decrease of iV" sh (r] c ) was found to be quicker in 
the strong-KNO-scaling-violating scenario than in the QCD-inspired one and more 
pronounced in larger intervals than in smaller ones. 

Finally the average number of clans increased with the width of the rapidity 
interval; the average number of particles per clan n CjS h also increased with the width 
of the rapidity interval and, contrary to N^, increased also with cm. energy. 

This situation raised intriguing questions, on when the asymptotic limit iV s h — > 1 
will be reached and on the consequences of its eventual earlier occurrence, say at 
LHC energy. 
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In (l| we concluded that an early occurrence of the limit iV s h — > 1 in FPS is of 
particular interest and leads in our framework to the onset of a new class of events 
with quite remarkable properties controlled by /c t h < 1, the effective benchmark of 
the new class of events. 

In particular are expected: 

a) very strong forward-backward (FB) multiplicity correlations (the FB multi- 
plicity correlations strength is close to its maximum, 6fb ~ 1) with leakage param- 
eter controlling the flow of particles from one hemisphere to the opposite one 0] 
p ~ 1/2, again its maximum value; 

b) an enhancement of two-particle correlations determined by the very large 
value of k^; 

c) a characteristic elbow structure in P* otal for large n and a narrow peak for n 
close to zero. Both trends are consequences of the log-convex gamma shape of the 
n-charged particles multiplicity distribution of the new component which shows a 
high peak at very low multiplicities and a very slow decrease for large ones, a general 
behaviour to be contrasted with Pythia Monte Carlo calculation predictions, which 
show at 14 TeV a second shoulder to be added to the first one already seen in the 
GeV energy domain. 

The situation is summarised in Fig. ^ within the framework of the weighted 
superposition mechanism of the three classes of events and in Fig. El for Pythia Monte 
Carlo predictions. It should be noticed that the total MD also in the case of Pythia 
can be fitted in terms of the superposition of three NB MD's with good chi-square 
per degree of freedom (107/66). The corresponding characteristic NB parameters of 
the two cases are given in Tabled andQ: respectively. Striking differences between 
the behaviours of events of the third class, seen just by inspection of Fig.sHJandEl 
are shown explicitely in the mentioned Tables. 



2 Perspectives: global properties in rapidity intervals 

Tevatron data seem to favour, among our scenarios, the one based on a strong 
KNO-scaling violation (i.e., with fc t ~J al — Ins). In fact it was shown by CDF p| 
that, in their data, the component rich in mini-jets violates KNO scaling in small 
rapidity intervals, although a cut-off at low pr has been used; in addition, E735 
data in FPS |6J] show again that the closest scenario is the one mentioned above, 
although discrepancies with previous UA5 results 0] are noticed in E735 results at 
lower cm. energies. Therefore, we decided to discuss in rapidity intervals the case 
of the strong-KNO-scaling-violating scenario only. 

In going from FPS to (pseudo)-rapidity intervals our main concern has been to 
be consistent with the scenarios explored in FPS. The weight of each component is 
obviously the same as in FPS. This leads to a particle density which shows an energy- 
independent plateau around rj = in both the soft and the semihard components, a 
constant k so f t value and a linearly increasing fctotai- It should be remarked that the 
plateau of the semihard component is less wide than that of the soft component, 
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Figure 1: Full phase-space multiplicity distribution for the scenario described in the text (solid 
line); the three components are also shown: soft (dashed line), semi-hard (dash-dotted line) and 
the third (dotted line). The inset shows a magnification of the low- multiplicity peak. 

but it is higher. See 0] for a detailed discussion on these points. 

For the third component, we have allowed for two extreme behaviours: (i) the 
third component is distributed uniformly over the whole of phase space and (ii) 
the third component has a very narrow and tall plateau and falls entirely within 
the interval \r]\ < 0.9. These two extremes are represented as a band in the figure. 
In the first case, the value of k t h has again been determined from the asymptotic 
behaviour of the average number of clans in the second (semihard) component, where 
it is a fixed fraction of the same quantity in FPS. Results are shown in Figure El and 
Table Cb. 

It should be pointed out that the general trend of P n vs. n is quite similar to that 
already seen in FPS: the elbow structure is clearly visible in both extreme cases; 
on the other hand, the narrow peak at very low multiplicities (also due to the third 
component) is hidden by the standard peaks of the soft and semi-hard components 
which are shifted to lower multiplicities than in FPS. 

This behaviour is clearly seen in Tables and by comparing n s h with n so f t 
and fc s h with k so f t in FPS and in rapidity intervals. In addition it turns out that 
iV soft (FPS) is slightly larger than iV sh (FPS) and N so{t (\r]\ < 0.9) is slightly smaller 
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Table 1: (a) Parameters for the extrapolated component MD's at 14 TeV in FPS; (b) same as 
(a), but in the pseudo-rapidity interval |ry| < 0.9, with the two extreme scenarios: evenly spread 
over the whole rapidity range (i) and concentrated in the small interval (ii); (c) Parameters of NB 
fits to Pythia predictions in FPS; (d) same as (c), but in the pseudo-rapidity interval |?7| < 0.9. 



a 


FPS 


% n k N n c 




soft 
semi-hard 
third 


41 40 7 13.3 3.0 
57 87 3.7 11.8 7.4 
2 460 0.1212 1 460 


b 


\r]\ < 0.9 


% n k N n c 




soft 
semi-hard 
third (i) 
third (ii) 


41 4.9 3.4 3.0 1.6 
57 14 2.0 4.2 3.4 
2 40 0.056 0.368 109 
2 460 0.1212 1 460 


c 


FPS 


% n k N n c 




first 
second 
third 


42 ± 4 52.8 ± 0.2 11 ± 1 19.3 ± 1.4 2.7 ± 0.2 
56 ± 4 123 ± 3 2.7 ±0.2 10.5 ± 0.7 11.7 ± 1.2 
« 2 468 ± 7 23 ± 3 70 ± 7 6.6 ± 0.6 


d 


\r]\ < 0.9 


% n k N n c 




first 
second 
third 


52 ± 7 5.1 ± 0.2 2.8 ± 0.3 2.9 ± 0.1 1.75 ± 0.1 
47 ±7 19.1 ± 1.5 1.9 ± 0.5 4.5 ± 0.7 4.2 ± 0.7 
w 1 87.7 ±3.7 11 ± 4 24 ± 5 3.6 ± 0.8 



Table 2: Forward-backward multiplicity correlation strength both in our scenarios and in Pythia 

Monte Carlo calculations. 





Pythia 
FPS \rj\ < 0.9 


Our scenarios 
FPS \rj\ < 0.9 


first / soft 
second/ semihard 
third 

total (weighted) 


0.39 0.28 
0.53 0.48 
0.91 0.80 

0.75 0.69 


0.41 0.25 
0.51 0.45 
0.9995 0.997 (i) 
0.9995 (ii) 
0.98 0.92 
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Figure 2: Three-component NBD fit to the full phase-space multiplicity distribution predicted by 
Pythia Monte Carlo calculations at 14 TeV; the three component NBD's arc also shown. 



than iVgh ( 1 77 1 < 0.9), but the average number of particles per clan is both in FPS 
and in \r)\ < 0.9 much larger in the semi-hard than in the soft component, a fact 
which confirms the general trend that semi-hard clans are larger than soft clans. 

As far as the third component is concerned, we have the mentioned extreme 
situations according to our assumptions: 

i) a situation in which the single clan is uniformly spread over the whole of phase 
space (in this case only 37% of the clan is contained within the pseudo-rapidity 
interval \r]\ < 0.9), fc t h is even much less than 1 (A^h ~ 0.056) and nth ~ 40. 

ii) a situation in which the single clan is fully contained in \r)\ < 0.9, i.e., its 
characteristic parameters are the same as those seen in FPS. It is quite clear that 
particle density in rapidity in this second case is much higher than in case (i). 

In Fig. 0] are shown the three components and their superposition predicted by 
Pythia at 14 TeV in the interval \t)\ < 0.9. We have taken the results from 
where Pythia version 6.210 [8] was run with default parameters using model 4 with 
a double Gaussian matter distribution. This set-up gives a reasonable description 
of lower energy MD data. The same events where then analysed both in FPS and 
in the rapidity interval \rj\ < 0.9. The MD's thus found have been fitted with the 
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Figure 3: Multiplicity distribution in \rj\ < 0.9 for the scenario described in the text (solid line); 
the three components are also shown: soft (dashed line), semi-hard (dash-dotted line) and the 

third (dotted line). 

weighted sum of three NBD's (it was not possible to use just two NBD's), which 
have been taken as the three components. The parameters of the NBD fits can be 
found in Table UJi (x 2 /NDF = 34/42). 

As in FPS, the behaviour of the tail is completely different from that predicted 
by our scenarios. It should be noticed that the average number of particle per clan 
is quite small in all components although the average number of clans increases 
from the first to the third component suggesting Poissonian MD structure in each 
component. 

Coming to FB multiplicity correlations, it is important to stress that our previous 
calculations [H of the correlation strength, &fb, for the two components, 

, _ a iJ^ + i 1 - a i)r^ + l^i 1 - - ni) 2 
FB ~~ W rfi 1 Z Z — ' { > 

a i\+bi + i 1 ~ "Oirfe + W 1 ~ «i)( n 2 - ^i) 2 

where, for the 2-th component, fii is its average multiplicity, Di its dispersion, hi its 
FB correlation strength and «j its weight, should be extended to accommodate the 
third one. Accordingly, we derived the following formula for the overall strength for 
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Figure 4: Three-component NBD fit to the multiplicity distribution in \<q\ < 0.9 predicted by 
Pythia Monte Carlo calculations at 14 TeV; the three component NBD's are also shown. 

an arbitrary number M of components: 

EM bjD} x sr^ M s^ M I- _ - \1 

i=i a i l + bi ' 2 2^i=l ^j>i a i a j\ U i U j) , x 

Applied to the three components case, Eq (jSJ) leads to the results shown in Tabled 
where the leakage parameter in rapidity intervals has been taken to be the same 
as in FPS In particular it is clear that the FB correlation strength for the 
soft component is quite larger in FPS than in \r]\ < 0.9 (&FB, S oft(FPS) = 0.41, 
&FB,soft(M < 0-9) = 0.25) and still remains larger in FPS than in \rj\ < 0.9 for the 
semi-hard component (6pB,sh(FPS) = 0.51, &FB, s h(M < 0.9) = 0.45), although it is 
always much larger both in FPS and in \rj\ < 0.9 for the semi-hard component than 
for the soft one. The third component tends to saturate in all cases its maximum 
value, which is 1. The total FB multiplicity correlation strength resulting from the 
composition of the contributions of all classes of events is larger in FPS (6fb = 0.98) 
than in \r)\ < 0.9 (&fb = 0.92) but closer to its asymptotic value. 

FB correlations in Pythia do not show remarkable differences with respect to 
our scenarios except in the third component and in the total MD where they are 
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Table 3: Energy density and corresponding parameters for our scenarios and for Pythia Monte 
Carlo. The volume V = itR 2 t has been computed with proton radius R « 1.1 fm and formation 









time t sa 


1 fm. 








a 


our scenarios 


soft 


semi-hard 


third (i) 


(ii) 


total (i) 


(ii) 




dn/dy 

\ ±^T J I Ivlc V 1 

e (GeV/fm 3 ) 


2.5 
0.41 


7 

1.6 


20 

tJUU 

4.7 


230 
54 


10.8 
2.5 


19.2 
4.5 


b 


Pythia 


first 


second 


third 




total 






dn/dy 
(E T ) (MeV) 
e (GeV/fm 3 ) 


2.5 
350 
0.41 


9.5 
500 
2.4 


44 
500 
10 




12.5 
500 
3.0 





considerably smaller; in particular the difference with the third component is striking 
(0.8 Pythia, ~ 1 our scenarios). Stronger FB correlations at hadron level suggest 
an extraordinary stronger colour exchange process at parton level in the last case. 
In addition, Bjorken formula [9J for the energy density, 



e = - 



3(E T ) dn 
2 V ~d~y 



(4) 



y=o 



where (Et) is the average transverse energy per particle, V the collision volume and 
dn/dy the particle density at mid-rapidity, has been applied in order to compare 
its predictions on pp collisions with those on nucleus-nucleus collisions. Parame- 
ters of the formula and results are shown in Table 0k: dn/dy values correspond to 
predictions for our scenarios already seen in Table ^ lacking general expectations 
for the average transverse energy (Et), we used for the soft component the value 
measured at ISR and, in a conservative way, the value measured by CDF for the 
other components (to be intended as a lower bound, which leads to lower bounds 
for the energy density as well). 
It should be noticed that: 

a) the energy density for the semi-hard component in our scenario at 14 TeV 
s 1.6, i.e., of the same order of magnitude of that found at AGS at 5.6 GeV in 
Cu collisions (e ~ 1.7). 

b) the energy density for the third component in the spread out scenario is 4.7, 
comparable with the value of e recently measured at RHIC in Au+Au collisions 

(e w 4.6). 

c) the energy density for the third component in the other extreme scenario 
(high concentration) is ~ 54 even larger, being dn/dy much larger, than the LHC 
expectations for central Pb+Pb collisions (e > 15). 

On the other side, Pythia prediction for its third component (see Table Eb) is 
contained within our extreme predictions, and it is of the same order of magnitude in 



is 
O- 

i.c 
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the other components. Notice that the full minimum bias sample is also intermediate 
between our scenarios results. 

Of course our calculation of e is only indicative and should be taken with cau- 
tion. Although the use of Bjorken formula for pp collisions as well as the choice of 
parameters is rather doubtful, we consider our results quite stimulating because they 
suggest the possibility that the same characteristic behaviour of many observables 
seen at RHIC energies in AA collisions could be reproduced at LHC in pp collisions. 

The implications of the existence of a third class of events in minimum bias pp 
collisions characterised by the presence of one or very few clans, i.e., by the cor- 
responding NB parameter fc t h < 1, has been examined in pseudo-rapidity intervals 
accessible at LHC (\r]\ < 0.9) paying attention to total charged particle multiplicity 
distributions, forward-backward multiplicity correlations and energy density. As- 
suming that the guesswork at the basis of our calculations will be confirmed, more 
than try to draw conclusions one should ask two intriguing questions for future 
experimental and theoretical work: 

1. what new physical phenomenon is hidden behind the squeezing of one or very 
few clans in the rapidity interval \rj\ < 0.9 with much stronger forward-backward 
multiplicity correlations than the soft and semi-hard components, and with energy 
density comparable with that of nucleus-nucleus at RHIC? 

2. what would be the QCD counterpart of this new phenomenon? would it be a 
signal of parton saturation? 
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